A local bias approach to the clustering of discrete density peaks 
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Maxima of the linear density field form a point process that can be used to understand the spatial 
distribution of virialized halos that collapsed from initially overdense regions. However, owing to 
the peak constraint, clustering statistics of discrete density peaks are difficult to evaluate. For this 
reason, local bias schemes have received considerably more attention in the literature thus far. fn 
this paper, we show that the 2-point correlation function of maxima of a homogeneous and isotropic 
Gaussian random field can be thought of, up to second order at least, as arising from a local 
bias expansion formulated in terms of rotationally invariant variables. This expansion relies on a 
unique smoothing scale, which is the Lagrangian radius of dark matter halos. The great advantage 
of this local bias approach is that it circumvents the difficult computation of joint probability 
distributions. We demonstrate that the bias factors associated with these rotational invariants can 
be computed using a peak-background split argument, in which the background perturbation shifts 
the corresponding probability distribution functions. Consequently, the bias factors are orthogonal 
polynomials averaged over those spatial locations that satisfy the peak constraint. In particular, 
asphericity in the peak profile contributes to the clustering at quadratic and higher order, with bias 
factors given by generalized Laguerre polynomials. We speculate that our approach remains valid 
at all orders, and that it can be extended to describe clustering statistics of any point process of a 
Gaussian random field. Our results will be very useful to model the clustering of discrete tracers 
with more realistic collapse prescriptions involving the tidal shear for instance. 



PACS numbers: 98.80.-k, 98.65.-r, 95.35. +d, 98.80.Es 



I. INTRODUCTION 

In the biasing scenario introduced by [l[ , virialized ha- 
los form out of initially overdense regions with a linear 
density (extrapolated to the redshift of interest) equal to 
5 C ~ 1.686. Since then, this picture has received con- 
siderable support from observational data. Even though 
dark matter halos are extended objects, they form a spa- 
tial point process as far as their clustering is concerned. 
However, this essential feature has remained elusive in 
most theoretical descriptions of halo clustering, which 
assume that halos are a Poisson sampling of a more fun- 
damental, continuous halo density field #h( x )- 

The peak formalism first proposed by @, H| in a cosmo- 
logical context is interesting because it is a well-behaved 
point process. In this approach, virialized halos are as- 
sociated with maxima of the initial density field. The 
displacement from their initial (Lagrangian) to final (Eu- 
lerian) position can be computed upon assuming phase 
space conservation Q . Clustering statistics of these dis- 
crete density peaks display many of the features present 
in measurements of halo clustering extracted from N- 
body simulations. In particular, discrete _density peaks 
exhibit a ^-dependent linear bias factor [5 
exclusion 0, Q , and a linear velocity bias 
these predictions have recently been tested in numerical 
simulations [HI El : peaks of the linear density field ap- 
pear to provide a good approximation to the formation 
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sites of dark matter halos with M > M*. 

However, despite recent progress towards the computa- 
tion of peak clustering statistics [|[ and a formulation of 
peak theory within the excursion set formalism (l2l . [l3j , 
discrete density peaks lack a clear connection with the 
more conventional local bias schemes fl4| . in which ha- 
los are approximated as a continuous field. Furthermore, 
while in the local bias model the computation of halo 
correlation functions is straightforward (though there are 
ambiguities regarding the filtering scale etc.), in the peak 
formalism calculations are particularly tedious owing to 
the peak constraint [H, 0, 0, flfl . In the most compre- 
hensive analysis thus far, ref. [4| succeeded in comput- 
ing the peak 2-point correlation ^ p k(^) up to second or- 
der, including the Zel'dovich displacement. They showed 
that some of the first- and second-order contributions 
could be obtained from a peak-background split formu- 
lated in terms of conditional mass functions. In con- 
trast to most analytic models of halo clustering, which 
assume that the (fc-independent) bias coefficients are the 
peak-background split biases, they derived this equiva- 
lence from first principles. However, they could not de- 
termine the physical origin of the other second-order con- 
tributions. Moreover, the peak constraint is clearly too 
simplistic to describe the clustering of low mass halos. 
In this mass range, one should consider more elaborated 
constraints involving the tidal shear etc. In this regards, 
it would be very desirable to find a simpler way of com- 
puting the correlation functions of generic point processes 
of a (Gaussian) random field. 

In this paper, we suggest a simple, physically moti- 
vated prescription based on the peak-background split 
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to compute the correlation functions of generic point 
processes driven by homogeneous and isotropic Gaus- 
sian random fields. We argue that clustering statistics 
of such point processes can be reduced to the evalua- 
tion of correlators of an effective continuous overdensity 
which, in the case of discrete peaks, is a function of the 
local (smoothed) mass density field and its derivatives. 
Our approach combines in a single coherent picture peak 
theory, peak-background split, local bias and the excur- 
sion set framework. For sake of clarity, we will focus on 
the 2-point correlation function of initial density peaks 
as computed in Q to explain the fundamentals of our 
approach. 

The paper is organized as follows. Sec. HI] furnishes a 
brief summary of clustering in peak theory. Sec. Mil is the 
central Section of the paper, where we present the con- 
nection between rotational invariants, peak-background 
split and a local peak bias prescription. Finally, SeclTvl 
discusses the implications of our findings. 



II. CORRELATION FUNCTIONS FOR 
DENSITY PEAKS 



where = v^ci/o^) is the characteristic radius of 
a peak (and not the interpeak distance). The three- 
dimensional Dirac distribution Soiv) ensures that all ex- 
trema are included. The factors of theta function 9h{\3), 
where A3 is the lowest eigenvalue of the shear tensor Qj , 
and the Dirac delta 8rAy — v') further restricts the set to 
density maxima of the desired significance v' . 

The (disconnected) Appoint correlations p^jp (or joint 
intensities) of density maxima are defined as the ensemble 
averages of products of n p ^(y, i? s ; x), 
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= (n p ] s .(v,R a ,xi) x • • • x 



n pk (i/, R s ,x N )\ • (3) 



For the Gaussian initial conditions considered here, mul- 
tivariate normal distribution are assumed to perform 
the ensemble average. In the particular case N = 1, 
(n p k(^, i? s , x)) = n v k{v,R s ) is the average, differential 
number density of peaks of height v identified on the fil- 
tering scale R s Q, 



We begin with a short recapitulation of the computa- 
tion of correlation functions in the peak formalism. Let 
8 S be the linear mass density field smoothed on scale 
R s with a spherically symmetric filter. For convenience, 
we work with the normalized variables f(x) = -^S s (x), 
r/i(x) = ±di5 s (x) and Qj{x) = ^d l d 3 5 s {x). Here, v is 
the peak height or significance, and 



n P k(^, R s ) 
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dkk 2{n+1) P s {k) . 



(1) 



are moments of the power spectrum P s (k) = (\S s (k)\ 2 ). 
A Gaussian filter is frequently adopted to ensure conver- 
gence of all the spectral moments <r„ , but one should bear 
in mind that the peak height associated with dark matter 
halos is always computed with a tophat filter (see [l3| for 
details). The first few spectral moments o~ n can be com- 
bined into a dimensionless spectral width 71 = o\l(a§(j\ 
that takes values between zero and unity. 71 reflects the 
range over which the smoothed power spectrum P s (k) is 
significant, i.e. 71 ~ 1 for a sharply peaked power spec- 
trum whereas 71 ~ for a power spectrum that covers a 
wide range of wavenumbers. 

Correlations of density maxima of S s can be evalu- 
ated using the Kac-Rice formula 0, [I?} ■ The trick is to 
Taylor-expand r]i(x) around the position x p k of a local 
maximum. As a result, the number density of (BBKS) 
peaks of height v' at position x in the smoothed density 
field 5 S can be expressed in terms of the field 5 S and its 
derivatives: 

33/2 

n pk (^,i? s ,x) = _|detC(x)|fc[r/(x)] 6 H [\ 3 (x)} (2) 
x 6 D [v(x) - v'\ , 
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In the last equality, 14 = (27r) 3 / 2 i??is the typical 3- 
dimensional extent of a density peak 12( . The functions 

Gn (7ij7i^) ar G defined in Appendix iBl In particular, 
the ratio G^/Gq 1 ^ is equal to the fcth moment u k of the 
peak curvature u. Similarly, the reduced 2-point correla- 
tion function for maxima of a given significance separated 
by a distance r = |r| = |x2 — xi| is 



£ pk (^, R s ,r) 



(2) 

Notice that, in p pk , we have ignored the shot-noise term 
«pk<5r)(x2 — xi) that arises from the self-pairs as it mat- 
ters only at zero-lag (in the peak power spectrum how- 
ever, this contributes a constant Poisson noise l/fi p k at 
all wavenumbers). 

The calculation of Eq.© at second order in the mass 
correlation and its derivatives is quite tedious 0, HI, HH 
because one must evaluate the joint probability distri- 
bution for the 10-dimensional vector of variables = 
(r]i(x a ), v(x a ), fyi(x a )) at two different spatial locations 
x Q = xi and X2, i.e. a total of 20 variables. Here, the 
components £4, A = 1, • • • , 6 symbolize the independent 
entries ij = 11,22,33,12,13,23 of Cij- Fortunately, as 
was shown in j4|, most of the terms nicely combine to- 
gether, so that the final result can be recast into the 
compact expression 
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The functions Q (r) are quantities analogous to a 2 ^ but 
defined for a finite separation r, 
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dkk 2(n+ Vp s (k) 3e (kr) , (7) 



where j^(x) are spherical Bessel functions. In the right- 
hand side of Eq.©, all the correlations depend on the 
filtering scale and the separation. However, the first line 
contains terms involving the first and second order peak 
bias parameters bi and bu (to be defined shortly), the 
second line retains a ^-dependence through the function 
1 + (2/5)d a lnG Q) ( 7l 

)7i J/ )|a=i solely, whereas the last 
two terms in the right-hand side depend on the separa- 
tion r (and R s ) only. Hence, unlike standard local bias 
expansions (Eulerian or Lagrangian), the peak 2-point 
correlation also exhibits quadratic terms linear in the 
second-order bias bn. These terms involve derivatives 
of the linear mass correlation £q and, therefore, vanish 
at zero lag. Clearly, they arise because the peak correla- 
tion also depends on the statistical properties of rji and 
Cv- 

Rcf . Q also showed that the Lagrangian peak bias fac- 
tors bjq(ki, . . . ,k^) can be constructed upon averaging 
over the peak curvature products of b v and b u , where 



b v (v, R s ) = 
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7i u 



b u (v,R s ) = — 

02 V 1 
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(8) 
(9) 



For peak of significance v on the smooting scale R s , the 
first order bias bi is defined as the Fourier space multipli- 
cation (we omit the dependence on R s and v for short- 
hand convenience) @ 



bi(fc) = &io + &oi£: 2 
where b\§ = b v and 



(10) 



&oi = K 



The overline designates the average over the peak curva- 
ture. 601 can be quite large for moderate peak heights. 
In the high peak limit v 3> 1 however, it is negligible so 
that bi(fc) is nearly scale- independent (like in local bias 
models). Similarly, the Fourier space expression of the 
second order peak bias bu is Q| 



bn(fci, fc 2 ) = 620 + bu(kl + kl) + 



b 02 k\k\ , 



(11) 



where k\ and k 2 are wavemodes and the fc-independent 
coefficients 620, 611 and 602 are 



b 20 (u,R s ) = bl- 



1 



bn{v, R s ) = b v b u 



^(1-7?) 
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(l-7i 2 ) 2 
7? 



(12) 
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and 



b 02 {v,R s ) = bl - 
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^(1-7?) 
u 2 — 2^\vu + j 2 u 2 
(l-7i 2 ) 2 



(14) 



(l-7i 2 ) 



By definition, bfj acts on the functions ^^(r) and 
!;i n2 \r) as follows: 



Air 4 



dk\ I dk 2 k 1 



2(m+l) , 2(n 2 + l) 



x K l {k ll k 2 )Ps(k 1 )P s {k 2 )3t l {k 1 r)j i2 {k 2 r) . (15) 

As pointed out by 0, the piece b 2 ^ 0) + (l/2)^° ) b 2 I ^° ) 
can be thought of as arising from the continuous, deter- 
ministic, local bias relation 

= kio<5 s (x) - b 01 V 2 5 s (x.) + ^b 20 5 2 ( x ) (16) 

- fo„<5 s (x)V 2 <5 s (x) + ife 02 [V 2 (5 s (x)] 2 , 

where the bias factors bij are peak-background split bias 
factors that follow from expanding the conditional peak 
number density in a series in the small background den- 
sity perturbation Si . This expansion is local in the sense 
that, except for the filtering, it involves quantities evalu- 
ated at x solely. However, an essential difference with 
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the widespread local bias model [14[ is the fact that, 
when computing the ensemble average (^pk( x i)<5pk( x 2)), 
we must ignore all powers of zero-lag moments (such 



(<^(xi)^(x 2 ))) to recover £ pk (r) 



since o 



as, e.g., ctq 

the latter does not exhibit such contributions (this 'no 
zero-lag requirement' also arises in the derivation of the 
'renormalized' bias parameters of (HI)- All the terms 
in Eq. (|16p are of course invariant under rotations since 
<5 p k(x) transforms as a scalar under rotations. Clearly 
however, this series expansion is not the most generic 
Lagrangian expansion we may conceive of (see, e.g., [l9j 
non nonlocal Lagrangian bias). 

Notwithstanding these results, [1] did not succeed in 
finding a physical interpretation of the other second-order 
terms in the right-hand side of Eq.©, even though it 
was pretty clear that they - at least partially - arise 
from coupling involving the components of the gradient 
r\i and the hessian Qj. 



rj^ 1 ^ can be obtained by projecting out the scalar and 
vector parts of the correlation of the Cartesian compo- 
nents rji with the projection operator P = e + <g>e + +e_ 

e_. The rule of thumb is that (r]{ fj^ ') = (rj^ rj2~ s 

where rja 



77^ (x a ), vanish unless s 



(r?i : 



(±iU±i) 
'12 



0. We find 



(19) 



Here and henceforth, the subscripts "1" and "2" will de- 
note variables evaluated at position X! and x 2 for short- 
hand convenience. Similarly, the symmetric tensor 
can be decomposed into its trace and traceless compo- 
nents. 



III. A INTUITIVE INTERPRETATION OF £ pk (r) 

In this Section, we propose an intuitive, physically mo- 
tivated explanation of Eq. © that is grounded in the 
peak-background split argument We begin with a 
brief introduction to the helicity basis, which was used 
in [if to compute probability distributions of the density 
field and its derivatives at two different spatial locations. 



A. Probability density in the helicity basis 

The 2-point correlation function of initial density peaks 
is the ensemble average of n p k(yi)?T.pk(y2) over the joint 
probability density Pa(yi) Y2\ r), where y a = y(x Q ) are 
the values of the field and its derivatives at position 
x Q . In what follows, n p k(y) will also designate Eq.@. 
Following 0, we can decompose the variables yj = 
(jj,:(x q ), v(x a ), Ca(x q )) that appear in the joint probabil- 
ity density P 2 (yi , y 2 ; r) in the helicity basis (e+,r,e_), 
where 



30 — eg 



r/r, 



y/2 



(17) 



Cij — g u Oij + Cij 



(20) 



3 Cj 



(T) 



The variables u 



-trC 



-C| and C (S) = C (0) 



the longitudinal and transverse hclicity-0 modes, Q, 



are 

(V) 

0, 



are the components of a transverse vector, £ 

(T) 

whereas Q- is a symmetric, traceless, transverse ten- 



sor, = C 

variables arc 



0. Explicit expressions for these 



C<*> EE Is^Clm = \ (3 



5 lm ) Ci, 



(T) _ 



1 J 



(21) 
(22) 

(23) 



and f, eg, and are orthonormal vectors in spherical 
coordinates (r, 9, <fi). The orthogonality relations between 
these vectors are e±-e± = rr = 1 and e + -e_ = e±-r = 0, 
where the inner product between two vectors u and v is 
defined as u • v = mVi = uiv % . Unless otherwise stated, 
an over line will denote complex conjugation throughout 
Sec. ESI 

In this reference frame, we decompose the first deriva- 
tives as 



rj = n^r + Tj [+L) e + + t/- 1 ^- . (18) 
W = rj-r and T]^ 1 ' = r] e± are the helicity-0 and 



where S a b, V^ c and are the scalar, vector and tensor 
projections operators (see, e.g., (20|). We have intro- 
duced factors of and \/2/3 in the decomposition 
Eq. (f2"0)) such that the zero-point moments of the helicity- 
0, -1 and -2 variables all equal 1/5 (see Eq. (|24|) below). 
The helicity-1 components of and their complex con- 



jugates are given by C^ 1 ) 
and C (±1) = C 



(V) 



e± 



y/3t±fi£ij, whereas C (±2) = 
5±C»j and C ( 



±2) 



C (T) e 



p3 



Here, 77 
-1 components 



y3/2 e l ±e±Cij are the t w0 independent hclicity-2 modes 
(polarizations) and their complex conjugates, respec- 
The correlation properties of 77^ and tivcly. Hereafter designating (^ s )(x Q ) as Ca\ the cor- 
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relation properties of these variables are the following: 



1 



5 4o 



(2) 



1 /l 



t(2) 
SO 



-J2) _ -- f 

2 t( 2 ) , 3 

7^ 2 + 3V 4 



18 (2) 

35 44 



(24) 



12 



(2: 



X t(2) 

5 4 ° 



(cr 2) c? 2) > 

and (C^sf ) = (Ci^C^' )- An the other correlations 
vanish. Note that the covariances are real despite the 
fact that the helicity-1 and -2 variables are complex. 

While the average peak number density only depends 
on the matrix M of covariances at the same location, 
the computation of the peak 2-point correlation func- 
tion £ p k(?") and higher-order clustering statistics from 
Eq. (|3|) generally involve covariances of the random fields 
at different locations. For £ p k(7~)j the covariance matrix 
C(r) = (yyt), where y = (yi,y2), is a 20-dimensional 
matrix that may be partitioned into four 10 x 10 block 
matrices: the zero-point contribution M in the top left 
and bottom right corners, and the cross-correlation ma- 
trix B(r) and its transpose in the bottom left and top 
right corners, respectively. Expressions for M and B(r) 
in the helicity basis can be found in the Appendix of 



latter is obtained upon setting Q n ^ = whenever £ ^ 
in the expression of B(r). As a consequence, C(r) retains 
the correlations {v\v 2 ), (^1^2), (uiu 2 ) and parts of the 

covariances (77^ ifc ) an d (Ci™^ C^)- However, Q did 
not provide a convincing explanation for their observa- 
tion. We shall do it now. 

Firstly, it is pretty clear that, since the peak 2-point 
correlation is invariant under rotations of the reference 
frame, it should be possible to express it in terms of ro- 
tational invariants constructed from the variables v : r\i 
and Qj . The peak significance v and the trace u are two 
obvious candidates, but they are not the only ones. The 
vector r\ of first derivatives and the traceless matrix £ 
yield two additional invariants, i.e. the square modulus 
rj ■ r] and the trace tr(£ 2 ). In the helicity basis, these 
invariants can be written 



(25) 



tr(C 2 



£(0)^(0) 

8=1,2 ' 



(26) 



B. Rotational invariants 

Translational and rotational invariance implies that 
n p k does not depend on spatial position, and that £pkM 
be a function of the distance r solely. In this regards, [4[ 
noted that, although the covariance matrix C(r) in the 
helicity basis (fT?)) does not depend on the direction f of 
the separation vector r, it is not equal to the angular av- 
erage covariance matrix C(r) = (l/47r) Jdilj. C(r). The 



The 3x3 symmetric matrix actually provides a third 
invariant with respect to rotations: the determinant det£. 
However, as we shall see below in the discussion of the 
peak-background split, because this determinant only en- 
ters the peak number density n p k(y) and not the 1-point 
multivariate normal distribution Pi (y) (the argument of 
the exponential is quadratic in the variables), it does not 
contribute directly to the peak bias. This suggests that 
we look at the covariances of f] 2 (x.) and C 2 ( x )j where 
these arc defined as 



?? 2 (x)^7,(x).r,(x) = -^ J^L J^ SlJ k u k 2j 5M6 s (k 2 y^ + ^ 

d 3 fci fd 3 k 2 



C*(x) 



C^(x) 



2a 2 



(2tt) 3 J (2tt) 3 



SiiS 



1 - ^SijSim ) ki l ki j k 2 ik 2m <5 s (ki)<5 s (k 2 )e 



(27) 

i(ki+k 2 )-x ^ ^8) 



where 6 S (k) are the Fourier modes of the smoothed den- 
sity field. As we shall see in Sec. MI C[ the variables 
3r/ 2 and 5£ 2 are distributed as chi-squared (x 2 ) variates 
with 3 and 5 degrees of freedom, respectively. Using ei- 
ther the Fourier space expression of j? 2 = ?7 2 (x Q ) (not 
to be confonded here with a cartesian component of the 
vector rf) or the fact that only components of identical 
helicity correlate, it is straightforward to compute the 
following correlators (for illustrative purposes, Appendix 



[A] furnishes details of the calculation of (riff] 2 )): 



2 
37f 



(^) 2 +2(^) 



= 1 



2 2 

(75 a 1 



SI SI 



(29) 
(30) 



Ignoring the zero-lag contributions, these terms corre- 
spond exactly (up to a sign factor) to some of those en- 
tering the second-order contribution of £ p k(?") in Eq.®, 
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with (£5; 1 ' /2 ^2o£i 1 ^ 2 ' 1 ) being proportional to {">]\v 2 ) hi 
particular. The computations of correlators involving 
£ 2 = C 2 (x a ) proceeds in a similar way although, in this 
case, it is much easier to sum the correlations among 
equal helicity components. For instance, 



Owing to rotational invariance, this probability density 
is a function of v, u, rj 2 and Q 2 solely (see, e.g., [2l| for a 



systematic analysis of distribution functions of homoge- 
neous and isotropic random fields). The quadratic form 
Qi(y) that appears in the exponential factor reads 



1 



2 (Cfcf) 2 



(31) 



2E((c^ +s) > 2 + (cr ) c^ ) > 2 )} 

s=l,2 ' 



where we used the fact that {C^ ] C2 s) > = (Ci^C^)- Aftcr 
some algebra, we find 



r( 2 h 2 



10 



(€ 



(2)\2 



v(^i 2) ) 2 



and 



= 1 

(Cl4) = i 



ha\a\ 



3(# /2) ) 2 + 2 



(3/2)n2 



(i)e(i) 



(32) 

(33) 
(34) 



The cross-correlations of n\ and (f with U2 in place of ^2 
are identical except for the superscript (n), which should 
be replaced by {n + 1). Again, all these terms can be 
found among the second-order contributions in the right- 
hand side of Eq. (|6]) . 

Therefore, the actual dependence of £ p k(»") on the in- 
variants r? 2 (x) and C 2 ( x )j whose covariances involve the 
correlation functions Q with I 7^ 0, is the fundamen- 
tal reason for C(r) being different from the angle average 
C(r). Those correlations, which arise upon expanding the 
joint probability density at second order, eventually all 
nicely combine together (and with terms proportional to 

(£0 ) 2 ano - (Co 2 ') 2 ) to yield the second-order correlators 

(v 2 i4), {ChD etc. 

The question then arises of the calculation of the co- 
efficients of these quadratic terms in the peak 2-point 
correlation function. We already know that the coeffi- 

(n) (n r ) 

cients multiplying products of the form Co Co are the 
quadratic peak-background split biases associated to the 
scalars v and u. Does this hold also for the coefficients 
multiplying (rj 2 r)%), (C 2 ^f) etc. ? 



C. Generalizing the peak background-split 



Qi(y) = 



271 vu 



2(1-7?) 



2 V 



5 2 
2 



(36) 



so that exp[— Qi(y)] retains factorization with respect to 
(u,u), rj 2 and C 2 - Furthermore, since r]( s > (with s = ±1) 
and (with s — ±1, ±2) are complex random variables 
with mean and variance 1/3 and 1/5 respectively and 
since (K (0) > = (^? (0) ) = and (r/°)C (0) ) = at the 
same spatial location, the quantities 3ry 2 (x) and 5£ 2 (x) 
are independent x 2 -distributed variables with 3 and 5 
degrees of freedom, respectively (similar conclusions can 
be drawn for the distribution of the components of the 
deformation tensor, see [H, H^)- Therefore, the 1-point 
probability density can also be written as 

Pi(y)d 10 y = u) dvdu 

x x 2 (3 )? 2 )d(3 ?/ 2 )x 2 (5C 2 )d(5C 2 ), (37) 

where M(y, u) is the bivariate normal 



and 



2ity/l^ 2 



exp 



v 2 + u 2 — 2^\vu 



2(1-7?) 



2 fe / 2 r(fc/2)' 



r.fc/2" 



l e -»/2 



(38) 



(39) 



is a ^-distribution with k degrees of freedom. Note that 
the distribution of £ 2 is coupled with the last rotational 
invariant det( (23|. However, it can be easily checked 
that, upon integrating over the (uniform) distribution of 
det£, we obtain the x 2 -distribution x 2 (5C 2 )- 

Rcf. [|[ discussed how the peak bias factors bij can 
be derived from a peak-background split. They argued 
that, while the fc-dependent piece b n o is related to the 
nth order derivative of the differential number density 
ripk, derivatives cannot produce the bias factors 601, fen 
etc. multiplying the fc-dependent terms. For this reason, 
they considered a second implementation of the peak- 
background split [24|, HH in which the dependence of the 
mass function on the overdensity of the background is 
derived explicitly. However, it is possible to write all the 
peak bias factors bij as derivatives of N(v, u) rather than 
?ipk- More precisely, the bij are the bivariate Hermite 
polynomials 



The probability distribution -Pi(y) that is needed to 
compute n p k is a multivariate Gaussian of covariance ma- 
trix M, 



Pi(y)d w y 



1 



(2 7 r) 5 |detM| 1 /2 



(35) 



Hij(v, u) = N{v, u)~ 



d_ 

dv 



d_ 

du 



M{v,u) (40) 



relative to the weight N(v, u), further averaged over the 
peak curvature u. Therefore, they are peak-background 
split biases in the sense that they can be derived from 
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— /d lu yn pk (y)^(^ U )P 1 (y) , (41) 



the transformation v — > v + e% and u — > u + £2, where 
Ej and £2 are long-wavelength background perturbations 
uncorrclatcd with the (small-scale) fields i/(x) and u(x). 
Ref. [3| did not express the peak-background split this 
way because they considered the effect of a background 
perturbation after the integration over the peak curva- 
ture. In terms of the rotational invariants introduced 
above, we can write the bu as 



where it is understood that P\ (y) takes the form Eq. ([37]) 
and d 10 y = dv du d(3rf) e?(5£ 2 ) . Factors of 1/cro an d 
I/02 are introduced because bias factors are ordinarily 
defined relative to the physical field 5 S (x) and its deriva- 
tives rather than the normalized variables. In practice, 
the integral is most easily performed upon transforming 
the 5 degrees of freedom attached to 5( 2 to the shape 
parameters v and w and the 3 Euler angles that describe 
the orientation of the principal axis frame (sec Appendix 
E]for details). 

In [23l | , it was noticed that, in the presence of non- 
Gaussianity, the 1-point probability density Pi (y) can be 
expanded in the set of orthogonal polynomials associated 
to the weight provided by -Pi(y) in the Gaussian limit. 
The same logic applies to the peak bias factors. Namely, 
the bij are drawn from the orthogonal polynomials as- 
sociated to Af(v,u), i.e. bivariatc Hermitc polynomials. 
Therefore, we expect that r/ 2 (x) and C 2 ( x ) also gener- 
ate bias parameters, and that these are drawn from the 
orthogonal polynomials associated with ^-distributions, 
i.e. generalized Laguerre polynomials. These are defined 
as 



z x d n 



i\ dx n 



(42) 



and are orthogonal over [0,oo[ with respect to the \ 2 ~ 
distribution with k = 2(a + 1) degrees of freedom. The 
orthogonality relation can be expressed as 

°dxx a e-*L£\x)L£Xx) = r ( n + « +1 ) Smn . 

n\ 

(43) 

The first generalized Laguerre polynomials are Lq (x) = 

1 and L^(x) = — x + a+ 1. 

Given the correlator Eq. (|29[) . the term proportional 
to l/erf in the right-hand side of Eq.([6|) indicates that 
the first-order bias parameter xio associated with the in- 
variant a 2 r] 2 (x) = (V6 s ) 2 (x) is xio = -3/(2cr 2 ). The 
aforementioned considerations suggest that we define 
the /cth-order bias factor as the Laguerre polynomial 
(— l) k L^ 2 \3r] 2 /2) averaged over all the possible peak 
configurations, i.e. 



v 2k Xko 



(-l) fc 



d w yn pk (y)L 



(1/2)/ 3?? 



Pi(y) 



(44) 



Taking into account the peak constraint, the first-order 
bias factors thus is 



1 



Xio 



<r 2 n pk 



rf lu yn pk (y) -77 2 -- Pi(y) 



3 

2a 2 ' 



(45) 
(46) 



which is precisely what we were expecting. Similarly, we 
define the bias parameter xofc associated to the invariant 
£ 2 (x) as the ensemble average of the Laguerre polyno- 
mial L^ 2 \5( 2 /2) orthogonal with respect to the weight 
xi(5C 2 ). Namely, 



rr 2 K 

02 XOk = — 

?1pk 



(3/2) (K 



d lu yn pk (y) J Lr / '( -5- My) (47) 



Note that, although we use a single symbol Xij to des- 
ignate the bias factors derived from the x 2 -distributions, 
the variables rj 1 and £ 2 , unlike v and u, are uncorrelated. 
The first-order bias factor thus is 



X01 



(y) ■ (48) 



To evaluate the integral, we first express the measure 
d(5£ 2 ) in terms of the ellipticity v and prolateness w, 



so that £ can be written as £ 



3v 2 



(see Ap- 



pendix [B] for details). A multiplicative factor of £ 2 will 
arise upon, e.g., taking the derivative of exp(— 5a£ 2 /2) 
with respect to a. In the notation of [4(, our fac- 
tor of £ 2 precisely corresponds to their derivative term 
— (2/5)9 Q lnGQ Q '(7i, 71^) evaluated at a = 1 (see Ap- 
pendix [BJ. Taking into account the factor of Gq (71, 71^) 
in the denominator, xoi can eventually be written 



X01 



2a 2 



1 



^9 a lnG^ ) (7i,7i^) 



(49) 



The physical interpretation of this result is straightfor- 
ward: C 2 ( x ) is a scalar that describes the asymmetry of 
the peak density profile. In the high peak limit, xoi — ^ 
reflecting the fact that the most prominent peaks are 
nearly spherical (see Fig. 9 of Q). 

The physical origin for the appearance of these orthog- 
onal polynomials can be found in the peak-background 
split. Long- wavelength background perturbations locally 
modulate the mean of the distributions A/"(*7 u), X 2 (3?7 2 ) 
and xi(5C 2 )- The resulting non-central distributions can 
then be expanded in the appropriate set of orthogonal 
polynomials. In practice, it is convenient to introduce a 
shift or translation operator T e to describe the action of 
a background perturbation on the distribution of rota- 
tional invariants. For the scalars v and u, we define the 
shift operator as 



f c = cxp (-£i<9„ - e 2 d u ) , 



(50) 



where t\ and £2 are small perturbations to the peak sig- 
nificance and the peak curvature, i.e. v — > v + t\ and 
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u — > u + 62 . The action of T e on the probability density 
JV(v, u) is to shift the (zero) mean of v and u by — ei 
and — £2, respectively (the reason for the minus sign is 
that Hermite polynomials include a factor of (— l) 1 " 1 "- 7 ). 
A straightforward calculation gives 

M{u,u)-^f e M{u,u) = m*-ji,"-t2) (51) 

A! (v, u) 
= f{e 1 ,e 2 )e ei ' Tob '' +e2 ' T2K , 

where b v and b u are given in Eqs.© and and f(ei, e 2 ) 
is the exponential factor in J\f(v, u) with the replacement 
v — > E\ and u — y 62 . The last expression is a generating 
function of bivariate Hermite polynomials. On expanding 
it in the small parameters ei and £2, 



(/(ei,e 2 )e £lCT ^ +£2CT2b "|pk) 



E^M (52) 

i,3=0 



we recover the bias factors bij once the results are aver- 
aged over all locations that satisfy the peak constraint. 
Note that the bias parameter b^2g defined in |l8| bears 
the same physical meaning as our 601: both represent 
the leading-order response of the tracer abundance to a 
uniform shift in the curvature of the density field. 

For the quadratic variables rj 2 and ( 2 , Eq.(|42p suggests 
that we express the shift operator in terms of both x 
and d x . The definition is somewhat cumbersome because 
we must take into account not only the ordering of x 
and d x , but also the factor of T(n + a + l)/n! in the 
orthogonality relation Eq. p5]) . A sensible definition of 
T e for the variable x = 3?y 2 and 5£ 2 is 



E 

3=0 

= r(k) 



j!r(|fc + i) 



I*(\/-%d x x) 



\d x x) 



a/2 



(53) 
(54) 



where I a (x) is a modified Besscl function of the first kind 
and the symbol :: of normal ordering is borrowed from 
quantum field theory. In the present discussion, the nor- 
mal ordering is defined as 



:{d x x) n : f(x) = d2(x n f(x)) 



(55) 



where f(x) is some test function. With this definition, 
the action of T e on a ^-distribution with k = 2(a + 1) 



degrees of freedom is 



(56) 



OO 



E 



(-§) 



(d xX y 



E 



-x/2 



-di 



(1) 



-x/2 



-x/2 



-(a) 



XfeWE r ^ fc+ ^ 



This is precisely the Laguerre series expansion of a non- 
central x 2_v ariate derived in [26| (see Appendix [Cj) . 
Therefore, 



where x fc 2 (x;A) 1S a non-central x 2 -distribution with k 
degrees of freedom and non-centrality parameter A = e. 
The latter is defined as the sum of squares A = J^ i=1 /if) 
where \ii are the means of the random variables. We 
can now read off the bias factors from the expansion of 
[Xfc] _1 ^cXfc m generalized Laguerre polynomials. For in- 
stance, 

'x; 2 (5C 2 ; 



£5 



X^.(5C 2 ) 



pk 



r(l) 



+ 3 



2j 

2 XOj 



(57) 



Note that the more common generating function 



(1-6)" 



cxp 



XX 



2 1 



$>"4 a) (z) (58) 



n=0 



appears to bear little connection with the non-central 
X 2 -distribution. 

To better understand the reason why the peak- 
background split generates a non-central x 2_ distribution, 
we note that, owing to the relation H2k(x) ~ Lk(x 2 ) be- 
tween Hermite and Laguerre polynomials, we could also 
have defined xio and xoi as second derivatives of Pi(y), 
where y is now the vector (v, u, t^Ca) of independent 
normal random variables such that rf = J2j=i W and 
C 2 = Xm=i Cl- A little algebra shows that 

Xl0 ^^- [d w y n pk (y)J2(-^-) 2 Pi(y), (59) 



Another way of writing this formula would be to absorb 
a factor of l/y/2 in the definition of ry^ 1 ) (and thus ex- 
plicitly deal with complex normal distributions). Analo- 
gously, we have 



X01 = — d w yn pk (y) 



5 

E 

A=l 



LA 

^80 



Pi(y)- (60) 
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This suggests that the effect of a background pertur- 
bation on rj 2 and ( 2 can also be thought of as shift- 
ing the components of the first derivatives according to 
Vi ~~ * Vi + £ 3ii and those of C according to (a — > (a + £5 a- 
The small perturbations e^i and e^a need not be the 
same for distinct i and/or A. However, owing to in- 
variancc under rotations, only the length of the vector 
Z)Li( e 30 2 = e 3 and J2 5 A=i( e 5A) 2 = £5 matter. This is 
the reason why the background perturbation effectively 
shifts the respective x 2 -distributions to non-central \ 2 ~ 
distributions, with non-centrality parameter A = e 3 and 
A = £5. Note that it should be possible to formulate this 
peak-background split with the conditional peak number 
density n p k(V, R a \§±, Ri) in a large-scale region of over- 
density Si, like in However, one should then consider 
two long-wavelength perturbations rjf and £ 2 i n order 
to describe the effect of the background perturbation on 
?y 2 (x) and C 2 ( x )j m addition to Si (which suffices to de- 
scribe the effect of the background wave for both v and 
u since these variables are correlated). 

To conclude, [13] also pointed out that, even though 
there is no functional relation n p k = J-(S) for dis- 
crete density peaks, it is nevertheless possible to define 
renormalizcd bias parameters as the expectation values 
c n ~ (J 7 ). However, he did not compute them explic- 
itly, nor specified what is T (though it is pretty clear 
that it is related to n p k). Here, we demonstrated explic- 



itly that each of the combinations (is, u), rj 2 and £ 2 of 
rotational invariants generates a set of orthogonal poly- 
nomials which, upon taking the ensemble average over 
all the possible peak configurations, yields a set of bias 
factors. Furthermore, we showed that these bias parame- 
ters can be constructed from a suitable application of the 
peak-background split to the probability densities char- 
acterizing the invariants. We will now demonstrate that 
we can interpret the peak 2-point correlation as arising 
from a functional relation of the form S p ^ = ^(Ss, . . .). 



D. A local bias approach to £ p k(?") 

First, let us make sure that Q obtained the correct 
expression for £ p k(f)- Adding all the second-order con- 
tributions induced by rj 2 , Q 2 and their cross-correlations 
with v and u, our result differs from theirs in that the 
last term in the right-hand side of Eq. © appears to miss 

a multiplicative factor of 1 + (2/^>)d a \aG^ t \^i, 7i^)| a =i. 
Checking the calculation of [|| , we found the missing mul- 
tiplicative factor in their Eqs. (A50) and (A51), in the 
form of g(r)tr(£ 2 ). This term was fortuitously omitted 
in their final expression of ^^(r). Consequently, the cor- 
rect answer is 



(o)r;2 c (oh 



;(l/2)r t (i/2h 



^(£rV£ 2 J ) ( 1 + ^lnG^i,^) 



_5_ 

3 
2o~\ 



(^Y + + v(^T ( 1 + ^lnGHT!,^) 



10, 



j( 2 h 2 



18 



(2h2 



(61) 



(^) 2 + 2(^) 



(i)\ 2 



3(4 3/2) ) 2 + m dr2 T\{l + ^lnG^( 7l , 7l .) 

r 



We note that this omission has an impact only on the 
small-scale (r < 20 /i _1 Mpc) peak correlation displayed 
in Fig. 2 of their paper. Their results concerning the peak- 
background split, the gravitational evolution of £ P k(?") or 
the scale-dependence of bias around the Baryon Acoustic 
Oscillation are unaffected. 

Even though we cannot write down a relation of the 
form n p k = Tx{S), the peak correlation function up to 
second-order can nonetheless be thought of as arising 
from a local bias expansion <5 p k = J-x(S s , . . . ), i.e. 

M x ) = &io*a(x) - VV 2 <5 s (x) + i& 20 5 2 (x) (62) 
- b n S s ( x )\7 2 5 s ( x ) + i6o 2 [V 2 <5 s (x)] 2 
+ xio (V<5 s ) 2 (x) + -xoi [SdidjSs - <%V 2 <5 s ] 2 x) , 



provided that we ignore all the contributions involving 
moments at zero lag. Nonlocality enters through the fil- 
tering solely (which is the reason why we still call it a 
local expansion) . As emphasized in Q , it is important to 
realize that 5 p k is not a count-in-cell quantity. Counts- 
in-cells can generally be constructed using the void gen- 
erating function, see e.g. [IH |2j|, but it is beyond the 
scope of this paper to compute moments of the peak fre- 
quency distribution function. Since all the bias factors 
are peak-background split biases obtained from a suit- 
able average of orthogonal polynomials, one could try to 
write down an expansion in terms of orthogonal poly- 
nomials in the variables S s , V 2 <5 S etc. such that all the 
contributions involving moments at zero lag cancel out. 
We will explore this possibility in future work. Note that 
this idea was put forward for the first time by (30j . who 
considered correlations of regions above threshold as a 
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proxy for luminous tracers. More recently, [3l| proposed 
an algorithm based on Hermite polynomials to extract 
/c-dcpcndcnt bias factors from cross-correlations between 
the halo and Hermite-transformed mass density field. 

To make connection with the formalism of [27| (see also 
[13), we define the Fourier space peak bias parameters 
c„(ki, . . . , k„) as the sum over all the contributions to 
<5 p k(x) from a given order. We thus have 



ncntial factor, we find 



ci(k)= {b w + b m k 2 )W(kR s 



(63) 



and 



c 2 (ki,k 2 ) = j&2o + &n (*? + *f) 



+ b m k\k\ (64) 



2xio (ki • k 2 ) + xoi 



3(k!-k 2 ) 



k^k 2 



W(k 1 R s )W{k 2 R s ) 



where W(kR s ) is the smoothing kernel. These definitions 
are consistent with those of the "renormalized" bias pa- 
rameters introduced by [13] who argued that, owing to 
rotational symmetry, the peak bias parameters should 
take the above functional form. In particular, the cor- 
respondence between the bias factors associated with rj 2 
and ( 2 and those of [27| is E 2 = 3x i, C 2 = ~^Xw an d 
D-2 = &02 — Xoi- We stress, however, that the peak bias 
factors discussed in this work have not been obtained 
by means of a renormalization procedure. We speculate 
that the local bias expansion Eq. ([62"T) can be extended 
to all orders to match the exact peak 2-point correlation 
function at all separations. Namely, 



°° 1 f d 3 ki 
x <5(k 1 )...(5(k„)e 4(kl 



c„(ki 



(65) 



where c n (ki , . . . , k„) is a sum over all the possible combi- 
nations of rotational invariants involving exactly n pow- 
ers of the linear density field S s and/or its derivatives. 

Our peak-background split approach provides a sim- 
ple way of predicting the bias coefficients associated with 
any rotational invariant quantity. The Hermite weight- 
ing scheme introduced by (3l| furnishes a practical way 
of measuring the biases bij from simulations. Clearly, 
their scheme could be extended to also measure the biases 
Xij- However, because discrete density maxima require 
a somewhat more sophisticated treatment of counts-in- 
cells, we leave this for future work. Here, we merely 
establish a recursion relation between the 6y by consid- 
ering either the property 



d_ 

dv 



■7i 



d_ 

du 



i(y) 



-ve 



-Oi(y) 



(66) 



or the generating function in Eq. (|51[) . In the latter case, 
upon substituting b u = {<Jq/ <Ji) 2 {v / <Jq — b v ) in the expo- 



5 c boi = ( ~ ) { v% _ dcbio) 



(67) 



at the first order whereas, at the second order, we obtain 



Sibu = Si b v b u 



-v + v S c b 1Q - S r b 2 o 



(68) 



and 



S 2 M = Sib 2 



^ Si 



^i(i-7i 2 ) 



(69) 



(Jo) (l-7i) 



v 2 {v 2 + 1) - 2u 2 5 c b w + S 2 b 



■211 



At least for n = 1,2, bjy with 1 < I < k < n can be ex- 
pressed as a linear combination of bko, 1 < k < n, plus a 
polynomial in v. This agrees with the fin ding s of [3l[ ob- 
tained within the excursion set approach [31j (except for 
the multiplicative factors of (ao/cri) 21 )- Similar relations 
should hold at any order in the peak bias parameters 
bij. This structure arises from the fact that the peak- 
background split acts on probability densities, which are 
continuous functions of space. For discrete tracers such 
as density maxima, the background perturbation affects 
the fields appearing in the Gaussian multivariate Pi (y) , 
but not those entering the expression of the peak number 
density n p k(y). The peak constraint weights the peak- 
background split series expansion such that the peak bias 
factors are recovered. 



E. Correlation functions for excursion set peaks 

To make connection with the clustering of dark mat- 
ter halos, we must ensure that the density in a tophat 
region centered at the peak location never reaches the 
colla pse threshold (S c ) on any smoothing scale R > R s . 
Ref. [33[ showed that enforcing the conditions S(R S ) > S c 
and S(R S + Ai? s ) < S c as in [34[ provides a very good ap- 
proximation to the first-crossing distribution when the 
stochastic walks generated from the variation of R s are 
strongly correlated. An important consequence of this 
result is the possibility of restricting the excursion set to 
those locations that meet the peak constraint [Hj]. The 
number density of dark matter halos per unit mass and 
volume is usually written 



n(M) 



JLf( v )— 



(70) 
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FIG. 1: First-order (left panel) and second-order (right panel) Lagrangian peak bias factors derived from the excursion set peaks 
mass function Eq. (|72[) assuming a moving barrier B(/jq) — 5 C + 0.43oo (see text). The bias parameters have been multiplied 
by the appropriate factors of o;, so that they all are dimensionless and multiply combinations of normalized variables in the 
effective local bias expansion Eq. (|62[) . In particular, the second-order bias factors 0-9620, 0002&11 and 0-2602 are associated with 
i/(x) and u(x), whereas afxio and o^Xoi weight the contributions from ?? 2 (x) and C 2 ( x )> respectively. Vertical lines mark the 
peak significance at which the halo mass is M = 10 12 , 10 13 , 10 14 and 10 15 M Q /h (from left to right). Dashed curves indicate 
negative values. Results are shown at z — for a ACDM cosmology with normalisation og = 0.82. 



where f(v) is the multiplicity function. Following the 
approach of [H, HH , the number density of peaks identi- 
fied on the filtering scale R s and satisfying the aforemen- 
tioned conditions is 

nESp(R s )dR s 

= ^rJ fe) j dW y n My)uPi{y)R s dR s . (71) 

where, for simplicity, we have assumed that the smooth- 
ing kernel is Gaussian, but it is straightforward to 
generalize these results to arbitrary filters. Therefore, 
dv/dR s = vR s (oi/<tq) 2 and we can write the excursion 
set peaks multiplicity function as 

/espM = (j) n ESP (R s )^- (72) 

\/2tt \V*J li u 
= Vn ESP (v) . 

Here, V = M/p\s the Lagrangian volume associated with 
the filter (usually tophat). /esp is the fundamental in- 
gredient in the mass function prediction of j 1 31 ] since it 
can be interpreted as a multiplicity function. It is pretty 
clear that, in the peak 2-point correlation, the constraint 
8{R S ) > 5 C and 8{R S + AR S ) < S c will translate into an 
extra multiplicative factor of (02 / <Ja)uR s dR a in the inte- 
grand. Therefore, the ESP correlation functions can also 



be obtained in the exact same way as that of the BBKS 
peaks, but with a peak number density 

n ESF (u',R s ,x) = ^ (-^j) |detC| Mr?) (73) 
Rf V 71^ / 

x e H {\ s )5 D {v-v') . 

The bias parameters for excursion set peaks are thus 
given by 

hi = — / d 10 y n ESP (y)i^ {u, u) P 1 (y) (74) 

cr l cr J 2 n E sp J 

x k0 = 4^- A 10 y-ESP(y)4 1/2 f^) My) 

<jf "-ESP J V 2 / 

x 0k = 4t^- A 10 y-ESP(y)4 3/2) (^) My) • 

"ESP J V 2 / 

They are similar to the bias factors of BBKS peaks, ex- 
cept for the fact that the fcth order moment u k of the peak 
curvature must be replaced by u k+1 ju in Eqs. (fT2"]) - (|14p . 
and Gq must be replaced by G± in Eq. (|49|) . For large 
values of lo = jif, G^ asymptotes to G^ w a~ 5 / 2 oj 4 . 

This implies that daluG^i^ji, to) converges towards -5/2 
in the limit uj — > 00. Hence, the second-order bias in- 
duced by the asymmetry of the peak profile also vanishes 
in the high peak limit for the ESP multiplicity function, 
as it should be. 
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To illustrate the behaviour of the ESP bias parameters, 
we follow the methodology adopted in Namely, we 
smooth the density field with a tophat filter, while stick- 
ing to the Gaussian filter to define rji and (so that 
the second spectral moment 02 remains finite). There- 
fore, we have 71 = <xf x / '(ctot^g), where the subscript 
T and G refer to tophat and Gaussian filtering, respec- 
tively, and x denotes a mixed filtering, i.e. one filter is 
Gaussian and the other is tophat. Next, we construct the 
mapping between and Rq by finding the Rq(Rt) for 
which (SgSt) = (§t)- Finally, to account for departures 
from the spherical collapse approximation, we consider a 
moving barrier of the (square- root) form B = S sc + /3ao, 
where (3 = 0.43 [35j | . For simplicity, we will ignore the 
scatter around B even though it is quite substantial in 
the range of v we are interested in. Fig[T] shows the 
first- and second-order peak bias factors at z — for 
a ACDM cosmology with <j% = 0.82. Note that, while 
M = 10 13 M Q /h translate into a significance of v sa 1.2, 
the actual height of v = 1.2 density peaks is B/cjq ps 1.6. 

These results can be generalized to arbitrary filtering 
of the mass density field. For non-Gaussian initial con- 
ditions, there are a couple of subtleties which will be 
discussed elsewhere 14511 . 



IV. DISCUSSION AND CONCLUSIONS 

We have shown that the 2-point correlation function 
£ p k(r) of discrete density peaks can be computed, up to 
second order at least, from an effective local bias expan- 
sion in continuous fields that are invariant under rota- 
tions of the coordinate frame. This local expansion is not 
a count- in-cell relation in the sense that <5 p k(x) is merely 
an effective overdensity that can be used to recover the 
true Cpk(^) from a trivial evaluation of (<5 p k(xi)<5 p k(x2)). 
One of the consequences is that there only is one phys- 
ically motivated smothing scale: the Lagrangian radius 
R s of the halos. Yet another important difference with 
the widespread local bias model is that one shall ignore 
all the contributions from zero-lag moments in order to 
obtain the correct £ p k(0- 

All the bias coefficients can be derived from a peak- 
background split argument in which the background per- 
turbation shifts the zero mean of the 1-point probability 
distribution functions of the rotationally invariant fields, 
unlike essentially all the other peak-background split for- 
mulations which consider a change in the number den- 
sity of the tracers. Consequently, it is possible to de- 
rive bias factors from a peak-background split argument 
even if the variables are integrated over. The resulting 
probability densities can then be expanded in orthogo- 
nal polynomial bases. For the normally distributed peak 
height z^(x) and curvature w(x), these are bivariate Hcr- 
mite polynomials whereas, for the chi-squared distributed 
r? 2 (x) and C 2 ( x ), these are generalized Laguerre polyno- 



mials. The peak bias factors are then obtained upon 
averaging the appropriate orthogonal polynomials over 
all the spatial locations that satisfy the peak constraint. 

We have demonstrated that our simple local expansion 
reproduces the 2-point peak correlation function £ p k(0 
computed at second order by [|[ after a tedious expansion 
of the joint probability density P2 (yi , y2 ; r) . We believe 
that it should remain valid at higher orders. Further- 
more, because discreteness enters the calculation only 
when averaging the orthogonal polynomials, we specu- 
late that this local bias expansion combined with the 
peak-background split approach presented here can be 
generalized to describe the clustering of any point pro- 
cess of a Gaussian random field. The great advantage of 
our approach is that it circumvents the computation of 
P2 (yi , y2 ; ?~ ) j an d requires only the evaluation of P±(y). 

Our approach can be easily generalized to more so- 
phisticated constraints involving, for instance, the tidal 
shear dj dj $ (x) , where $(x) is the gravitational potential. 
As noted in [l^, [HJ , the quadratic invariant s 2 (x) (the 
equivalent of our £ 2 ( x ) but with 5 replaced by $) follows 
a ^-distribution with 5 degrees of freedom. Therefore, 
we expect that its associated bias parameters are given 
by some suitable average of the Laguerre polynomials 
L k If the (nonspherical) collapse occurs at the 

spatial location of density peaks or includes the depen- 
dence on the large scale environment, then the x\ dis- 
tribution will be replaced by the appropriate conditional 
probability density (36l . |37| , to which we shall apply the 
peak-background split in order to read off the new bias 
parameters. 

Corrections induced by nonlinear gravitational evolu- 
tion can also be decomposed into rotational invariants 
[38l - l40| . Therefore, if one ignores the diffusion kernels 
(i.e. the propagators introduced by [HI), then it is 
straightforward to find explicit expressions for the Eu- 
lerian bias parameters in terms of local and nonlocal La- 
grangian bias factors [ljl l42T - |44j . This procedure can 
clearly be applied to our effective bias expansion Eq. (152^) , 
with the important caveat that discrete density peaks 
exhibit a statistical velocity bias. Q. Notwithstand- 
ing this, we expect from the structure of the F2 kernel 
that the bias factors Xij remain constant with time, in 
agreement with the findings of Q|. For a more realistic 
treatment of gravitational motions, it should be possible 
to compute the evolved 2-point peak correlation ^ p k(^, z) 
in the framework of the integrated perturbation theory 
proposed by [27j. We leave all this to future work. 
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Appendix A: Computing correlators 

For illustation, we evaluate the cross-covariance (rjlrfe) 
of the square modulus of the gradient ?7 2 (x) at two dif- 
ferent spatial locations Xi and X2. We have 



(^ 2 )=^{n/0}( k 1 - k 2 )( k 3-k 4 ) 



(Al) 



x (5 s (k 1 )5 s (k 2 )(5 s (k3)5 s (k 4 ))e 



i(ki+k 2 )-X2+i( k 3+k4)-xi 



i + 4 



d3kl fd " k2 (k^) 2 ^)^) 



(2tt) 3 J (2tt) 3 



x e i ( k i+ k 2) r 

1 + (~sj fiijfiimJti(r)Jj m (r) , 



where 



Jij(r) 



d 3 k 



ki kj P s {k^)e 



(A2) 



(27r) 3 '" 1 r 

To evaluate Jij(r), we express ki in terms of the compo- 
nents ki = ki/k of the unit vector, and take advantage of 
the fact that the integral over the angular variables is 



1 

4tt 



jo(kr) + ji{kr) Sij - j 2 (kr)nfj . (A3) 



Therefore, the product SijSi m Ju(r)Jj m (r) becomes (we 
omit the r-dependence for conciseness) 



1 



:(!) 



(1) 



(^) 2 + 2(^) 



3 

(i)\ 2 



(A4) 
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which yields Eg . ([29)) once the multiplicative factor of 
(2/tTf) and the additive zero-lag contribution are ac- 
counted for. 



of asphericity expected for a peak, 



F(u, V, w) = (it - 2w) (u + w) - 9v 2 v (v 2 - iu 2 )(Bl) 



Appendix B: Shape factor for peaks 

In terms of the ordered eigenvalues Ai > A2 > A3 of 
the hessian matrix Qj, the asymmetry parameters that 
quantify the departure from a spherically symmetric peak 
density profile are v = (Ai — Aa)/2 and w = (Ai — 2A2 + 
A3)/2. The peak constraint together, with our choice 
of ordering, impose the four conditions v > 0, —v < 
w < v, (u + w) > 3v and u > 0. Following Q, we also 
introduce an auxiliary function that measures the degree 



This function scales as 



in the limit m> 1. 



In Q , the peak 2-point correlation up to second order 
(i.e. terms quadratic in the correlation of the density 
field and its derivatives) is written as the sum of the 

(2i) 

linear contribution and three second-order terms Q k , 

(22) 

i = 1, 2, 3. In particular, Q k contains all the terms for 
which the independence cannot be expressed as a poly- 
nomial in the linear and quadratic bias parameters bij. 
Their expression is phrased in terms of 




and its integral over the nth power of the peak curvature 
u times the w-dependent part of the one-point probability 
distribution, 

roo ( a; _ tl , ) 2/2(l_ 7l 2 ) 

Gi a \ Ju w)= dxx n f(x,a) . (B3) 

Jo V 27r ( 1_ 7f) 

These functions are very similar, albeit more general than 
those defined in Eqs (A15) and (A19) of Q. 

Appendix C: Non-central chi-squared distributions 

The probability density of a non-central chi-squared 
distribution x^ 2 (x; A) with k degrees of freedom and non- 
centrality parameter A is given by 

p -{x+X)/2 /tX q/2 / \ 

x k 2 (x; A) = — - — (j) i a (Va^J , (Ci) 



where a = k/2 — 1 and I a (x) is a modified Bessel func- 
tion of the first kind. Ref. [26| proposed the following 
Laguerre polynomial expansion, 



j=0 \ 2 J I 

The non-central x 2 -distribution can also be represented 
as a Poisson-weighted mixture of central x 2 -distributions 
(this was used by, e.g., [n| to estimate the nonlocal La- 
grangian bias induced by ellipsoidal collapse). Note, how- 
ever, that this representation does not make apparent the 
connection with the bias parameters Xij ■ 



